Group velocity of the acoustic eigen-modes in sonic crystals 
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In this brief report, the group velocity of the eigen-modes in sonic crystals is derived, and shown 
to equal the averaged energy velocity of the eigen-modes. How the group velocity can be used to 
describe acoustic energy flows in sonic crystals is discussed. 
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The phenomenon of band structures, revealed as waves 
propagate through periodic structures, was put on a solid 
foundation in the context of Bloch's theorem^]. It was 
first studied for electronic systems. Since the central 
physics behind electronic band structures lies in the wave 
nature of electrons, it is natural to extend the study 
to any wave systems in general, and to electromagnetic 
(EM) and acoustic systems in particular. 

Not excluding other efforts, band structures were ear- 
lier addressed for general waves in periodic structures by 
Brillouin0, then to EM waves by Yariv and Yehj3|, later 
by YablonovitchQ and John[j|. Since then, the study of 
EM waves in periodic structures has been booming, even- 
tually extended to the exploration of acoustic waves in 
periodic structures (e. g. Refs. 0, 0, Hi HI ) i leading to the 
establishment of the field of photonic crystals (PCs) and 
sonic crystals (SCs). 

Sonic crystals (SCs) are made of periodically struc- 
tured materials which are sensitive to acoustic waves, and 
have been studied both intensively and extensively. The 
exciting phenomenon of band structures in sonic crystals 
allows for many possible applications. It has been recog- 
nized that SCs could be used as sound shields, acoustic 
filters, acoustic flow guiding, sonic crystal lenses, and so 

One of the top issues in the research of acoustic crys- 
tals is how to describe the propagation of acoustic waves 
in crystal structures, by analogy with the photonic crys- 
tals. The common theoretical approach to electromag- 
netic propagation in periodic media has been given in 
Ref. 0], and may be summarized as follows 16] . The 
Maxwell equations are first derived for waves in periodic 
media. By Bloch theorem, the solution can be expanded 
in terms of Bloch waves. The solution is then substituted 
into the governing equations to obtain an eigen-equation 
that determines the dispersion relations between the fre- 
quency and the wave vector that lies within the first 
Brillouin zone for the eigen-modes. These relations are 
termed as frequency band structures. Since it has been 
proved 0,0] that the averaged energy velocity v e equals 
the group velocity which can be obtained as the gradi- 
ent of the dispersion relations with the respect to the 



space of wave vectors, i. e. v g = V^u, the investigation 
of electromagnetic propagation in periodic structures is 
thus reduced to the calculation of the group velocity from 
the band structures. 

Simply due to similarities between the electromag- 
netic and acoustic waves, it is naturally expected that 
the above prescription can be automatically adopted for 
sonic crystals, and this recipe has been indeed used for 
SCs in the literature (e. g. Ref. 0])- To the best of 
our knowledge, however, so far there has been no at- 
tempt in the literature to prove this analogy and con- 
nection between PCs and SCs. Here, we wish to bridge 
the gap. In this paper, we will provide a proof that in 
analogy with PCs, the acoustic energy velocity of the 
eigen-modes of a sonic crystal averaged within a unit cell 
also equals the group velocity of these eigen-modes, de- 
fined as v g = V^w. Moreover, we will discuss how the 
group velocity can be used to describe acoustic energy 
flows in sonic crystals. In addition, following the work on 
photonic crystals [l9|. it will be shown that measuring the 
unit-cell volume averaged acoustic energy velocity can be 
reduced to the measurement of the corresponding unit- 
cell surface averaged components. Such a reduction is ex- 
pected to facilitate the practical measurement of acoustic 
energy flows. 

General - In the linear region, the propagation of 
acoustic waves are governed by three physical princi- 
ples: Newton's second law, conservation law, and equa- 
tion of states. And acoustic waves are characterized by 
three major physical quantities: variations in pressure 
and mass density, and the displacement velocity, denoted 
respectively by 8p, 8p, and u. The governing equations 
are then 



_ dSp du 

These equations lead to 
V • u ■■ 



--V(Sp), and (Sp) = i(djp). 
p & 



1 dSp 
c 2 p dt 



So we have two main equations 
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V-u = 



1 dp 
c 2 p dt 



(1) 
(2) 

(3) 
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and 

du 1 

8i = - P VP > (4) 

where we ignore i S\ 

From Eqs. J2J| and (QJ, we have the wave equation 



V • ( -Vp 



1 <9 2 p 
c 2 p d 2 t 



0. 



(5) 



Multiplying Eq. © with p*, and multiplying Eq. Q 
with it*, and then adding the correspondingly the com- 
plex conjugated equations, we will get 

-V.(up* + tf* P ) = -^ -^-f^f (6) 



Here we assume that the band structure is known. Read- 
ers may refer to Ref. [jj for details. We will prove that 
the acoustic energy velocity of the eigen-modes of a sonic 
crystal averaged within a unit cell equals the group ve- 
locity of these eigen-modes. 

Considering monochromatic waves in sonic crystals, 
from Eqs. © and l@J 



V • U = 



ILOpiL. 



Vp 



(8) 
(9) 



Hereafter we drop out the time factor e %u}t . 

According to Bloch's theorem, waves of the eigen- 
modes in periodic structures can be expressed in terms 
of Bloch functions, 



This equation can be written as 



d 



(7) 



iK- 



AK- 



P = P K e "~ ' ■ ( 10 ) 
Taking these two equations into Eqs. (JSJ and 10, we have 



with J = Re[p*u], and U = 2^\p\ 2 + f \u\ 2 - Eq. 

represents the energy conservation. The quantities J and 
U are the acoustic intensity and energy respectively. 

Proof of V guj = v e - Now we consider sonic crystals. 
We will consider three dimensional cases; the extension 
to two dimensions is straightforward. The band struc- 
tures can be readily computed by the plane- wave method. 



luJ 



K 



Vpg + iKpj 



K 



c*p 

ILUpUg. 



(11) 

(12) 



Considering an arbitrary infinitesimal change in K, 
which induces the changes in u, u, and p, Eqs. (flip and 
(fT2"|) become 



J 



V • Sug + iSK ■ ug + iK • Sug 



—Sivpg 
c £ p 



-z-uSpg, 
c z p 



VSpg + iSKpg + iKSpg = ipSujug + ipuSug. 



(13) 
(14) 



The corresponding complex conjugates are 



V • 5u*g - iSK ■ u* g 



—* % % 

iK ■ 5ug = -—Sujpg - —ujSpg, 



WSp*g - iSKp R - iKSp*g 



c^p " c p 
-ip6u>u*~ — ipuSu*- 



A" 



(15) 
(16) 



Multiplying Eq. l [T5jl with p%, and Eq. l(l5|) with pg, equations, we have 
then making a subtraction between the two resulting 



2ilm[p*gV ■ Sug] + 2i5K ■ Jg + 2iK ■ Re\p*g6ug] 



2% 
c z p 



2i 

5uj \p k \ 2 + ^MPR^Pg] 

o p 



(17) 
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where 



Jg = Re\p*gtig] 



Multiplying Eq. {H} with tig, and Eq. JT§J with tig, 



then making a subtraction between the two resulting 
equations, we have 



2ilra[tig ■ VSpg] + 2iSK ■ Jg + 2iK ■ Re{ugSpg] = 2ipSLu\ug\ 2 + 2ipujRe[v*g ■ Sug] 

I 



(18) 



Adding Eqs. (|17|l and (JTSJ , and taking into account Eqs. <|11|) and (|12|1 . we have 



2ipSu\ug\ 2 + 2iIm[V • {p*gStig) + V • (tigSpg) 



(19) 



r 



Due to the periodicity, the unit cell integration for any 
periodic function A can be shown to be zero, i. e. 
J c drV ■ A = 0. Now performing the unit cell integra- 
tion for Eq. I|19|l . we have 

(^bi?| 2 + fK'| 2 )^-^-(4) ; (20) 

where (^4) = y- J c dfA, with V c being the volume of the 
unit cell. 

We define the unit-cell volume averaged energy velocity 

as 



Then we have 



Since 



Suj = v £; • SK. 

e,i\ 



Suj = V gLU ■ SK, 



(21) 



(22) 



(23) 



we have finally proved the identity for the group velocity 



defined as v g g = V ' gto: 



(24) 



When the index K for the eigen-mode is ignored, we 
simply have 

v g = v e , (25) 

with v g = V guj, which was first derived for PCs. 

From the above, we know that V^w may describe the 
averaged energy flow in sonic crystals. For example, 



the direction of averaged energy follows that of V^a>. 
As such, in the study of photonic crystals, Vj^w is re- 
garded as the key quantity in discerning the photonic 
flows. Here, however, we must stress that the average 
of the energy velocity is performed over the whole unit 
cell of periodic media, under what condition such aver- 
aged energy flow can depict the actual energy flow in 
periodic media is worth considering. Moreover, since in 
actual measurements it is often hard to detect the phys- 
ical quantities inside the periodic media, how to obtain 
the group velocity without having to put a detector into 
the media needs also to be considered. In the following, 
we will consider these two questions. 

First, due to the periodicity, the local acoustic current 
of the eigen-modes can be written as 



G 



Jg{G)e 



iG-r 



(26) 



where G denotes the reciprocal vectors. It is shown that 

G 

= E4(g)^,g = 4(°)- ( 27 ) 



From this result, we clearly see that when and only when 
the component with G — dominates, the group veloc- 
ity v g = Vjjo; can represent well the actual wave flows. 
In this case, the eigen-modes will behave effectively as a 
plane wave in a free space. Some recent studies indicate 
that this may not be generally valid, and have caused 
in certain circumstances some confusions about whether 
waves are diffracted or refracted, referring to the discus- 
sions in Refs. [20I l2l| . Additionally, when a crystal is 
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ensonified by an external wave, many eigen-modes may 
be excited for a given frequency, adding complications to 
the determination of the actual acoustic propagation. 

Reduction of integration - Second, it can be shown 
that the unit volume averaged acoustic intensity (Jg) = 

Y~ J c drJg can be represented by surface averaged acous- 
tic intensity. This is important, since it will allow us 
to measure certain physical quantities at surfaces rather 
than putting detectors into the volume of samples. To 
briefly show this, we will take the procedure stemming 
from the excellent note [n|. For details, readers should 
refer to Ref. 0. 

We consider an acoustic lattice with the base lattice 
vectors a\, a^, and 03. The base vectors for the reciprocal 
lattice &i, 62 and 63 are determined by d{ ■ bj — 2n8ij. We 
will show that 

(■4) ■ k = 4" / dsjl i< ■ k (28) 

°i J Si 

where bi is the unit vector of bi, given by bi — b% , and 
Si denotes the three surfaces of a unit cell, characterized 
by S2 x CI3, (J3 x ai, ai x 0,2 respectively. 

For the stationary case, i. e. V • Jg — 0, we have from 
Eq. ((23) 

3-Jg{8)=0 (29) 

To show the identity in Eq. lj2"S|l . we consider the sur- 
face 5*3 as the example. Writing G — J2i=i m ibi, we have 




= J2 4(^363) -63. ( 3 °) 
Applying Eq. to Eq. (J20J), we arrive at 

±- f dsJg. ■ b 3 = 4(0) • 63. (31) 

A 3 JS 3 

From Eq. (|27|l and by comparing to Eq. I|31|) . we get 

(Jit) ■ bs = ^-J dsJg ■ 63. (32) 

Following the same procedure, we can have the similar 
results for the other two surfaces. Thereby we prove 

(4)'^ = tt / dsfg-bi. (33) 

Oi J Si 

This result indicates that the components of the volume 
averaged acoustic intensity can be reduced to the surface 
averaged counterparts. Clearly, this is expected to be 
very useful to actual measurements. 

In summary, we have provided a proof that in analogy 
with PCs, the acoustic energy velocity of the eigen-modcs 
of a sonic crystal averaged within a unit cell also equals 
the group velocity of these eigen-modes, defined as v g = 
Vjjw. We have also discussed the conditions that the 
group velocity can be used to describe acoustic energy 
flows in sonic crystals. Then it is shown that the volume 
averaged acoustic energy velocity can be reduced to the 
surface averaged components. 
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